BASICS OF QUANTUM MECHANICS, GEOMETRIZATION AND 
SOME APPLICATIONS TO QUANTUM INFORMATION 

JESUS CLEMENTE-GALLARDO AND GIUSEPPE MARMO 

Abstract. In this paper wc present a survey of the use of differential geo- 
metric formalisms to describe Quantum Mechanics. We analyze Schrodingcr 
framework from this perspective and provide a description of the Weyl-Wigner 
construction. Finally, after reviewing the basics of the geometric formulation 
of quantum mechanics, we apply the methods presented to the most interest- 
ing cases of finite dimensional Hilbcrt spaces: those of two, three and four level 
systems (one qubit, one qutrit and two qubit systems). As a more practical 
application, wc discuss the advantages that the geometric formulation of quan- 
tum mechanics can provide us with in the study of situations as the functional 
independence of entanglement witnesses. 



1. Introduction 

Interference phenomena of material particles (matter waves) as electrons, neu- 
trons, atoms, etc provide the most striking evidence for the need to elaborate a new 
mechanics which goes beyond and encompasses classical mechanics. 

At the same time, phenomena like photoelectric and Compton effects show a 
"corpuscular" behavior of radiation requiring, therefore, the need for a revision of 
the classical description of radiation. 

The associated quantitative results suggest that wave-like and corpuscular-like 
attributes satisfy the following Einstein- de Broglie relation 

p a dx a - Edt = h(k a dx a - uidt) (1) 

where (x a ,p a ) are the coordinates of the phase space and t stands for the time. This 
relation between the Poincare one form on the phase-space over the space-time and 
the optical phase-space establishes a relation between momentum (p a ) and energy 
(E) of the corpuscular behavior and the wave-number (k a ) and frequency (w) of the 
wave behavior. The dimensional proportionality coefficient is the Planck constant. 

We may use this relation to predict under which experimental conditions light 
will behave like a "corpuscule" and an electron will behave like a "wave". 

By means of this relation it is possible to conceive of interference experiments like 
the one of the double slit. When carried on with electrons, wc find some peculiar 
results for which we do not have a simple interpretation in the classical framework. 
The actual experiment has been performed in such a way that at each time only one 
electron is present between the source and the screen. It is found that the electron 
impinges on the screen at "given points" and leave a spot. 

All the spots are alike, there are no "half-spots" . After a few hundred electrons 
have passed we find a picture of spots erratically distributed on the screen, as we 
can see in the following picture 
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However, with several thousands electrons, we get a very clear interference pat- 
tern. 




Experiments of this kind have been performed (see [37l [T4j l28l 139] ). 

The same findings occur if we experiment with photons (light quanta), again 
with an experimental arrangement that makes sure that only one photon is present 
between source and screen at each time. 

These interference aspects call for a theory with a wave-like description, having 
a statistical-probabilistic character, along with an intrinsically discrete aspect (i.e 
no half-spots should be possible). 

We could try now to perform an experiment to watch which path the electrons 
follow on its way from the source to the screen. An experimental setup uses a 
conducting plate and the mirror image of the charge to find out which region will 
be "heated" by the passage of the electron (see [36]). The setup is as follows: 




We find that for z ss the interference pattern disappears: 
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We may try to summarize the main findings of our experiments: 

• electrons arrive at the screen in identical lumps 

• there is a wave associated with the electrons to describe the interference 
pattern 

• there is a probability for the electron to impinge on the screen within a 
preassigned region 

• under some appropriate experimental conditions we recover a classical-like 
behavior of the electrons, i.e. the interference pattern disappears. More 
specifically, the more visible the interference pattern is, the less distin- 
guishable the paths are. 

What we have described is quite counter-intuitive for particles, but it is even 
more startling for light. We have to accept that it is not so simple to elaborate a 
model capable of accounting for these results within the classical framework. Indeed 
light provides the most relevant example of the dual nature of quantum objects: its 
oscillatory properties were suitable to verify the electromagnetic theory of Maxwell, 
while its lumpincss (the photons) signalled the birth of modern quantum theory. 
The dualism between the two pictures may appear to be in contradiction, however 
it constitutes the first example of what is known as complementarity in quantum 
theory, i.e. the possibility to display both wave and corpuscular properties. 

From the historical point of view, things went differently because inconsistencies 
arose already in the derivation of the spectral distribution of energy density of a 
black body radiation. Max Planck proposed his well known ad-hoc modification of 
the standard phenomenological approach. His result, however, did not explain the 
physics behind the successful result. 

Another deficiency of the classical theory arose at the level of the models available 
for atoms and molecules, where it was not possible to account for the stability of 
atoms and molecules while explaining the detected atomic spectra. An account of 
the experimental background material can be found in most textbooks of Quantum 
Mechanics, for instance in ([17). 

Various efforts of the youngest theoretical physicists of that period gave rise 
to two alternative, but equivalent, formulations of quantum mechanics, usually 
associated to the names of Schrodinger and Heisenberg. As it is usually the case, 
though, the final theory grew out of the efforts of a reasonable large community of 
theoretical physicists during the first half of the XXth century. 

This paper is an expanded version of some lectures delivered by one of us (G. 
M.) at the "Advanced Winter School on the Mathematical Foundation of Quantum 
Control and Quantum Information" which took place at Castro Urdiales (Spain), 
February 11-15, 2008. The audience was composed of PhD students in Mathematics 
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and in Physics, and this explains the pedagogical nature of this paper. The content 
of the lectures has been revised and expanded in a substantial way by both of us. 

2. SCHRODINGER'S WAVE MECHANICS 

2.1. Introduction. Of course, if there are waves to describe the wave-like behavior, 
there ought to be a wave equation to describe the evolution of the system in time. 
This equation is called the Schrodinger equation and has the following structure: 

I*— s**- (2) 

The solutions, as functions of time, describe the time evolution of the initial condi- 
tions of the system. The complex valued function ip is called the wave function 
and is defined on the configuration space of the system we are considering. From 
a physical point of view it is interpreted as a probability amplitude, but this inter- 
pretation requires, for consistency, that 

ip*ipd/J, = l, (3) 

D 

i.e. the probabilistic interpretation requires ip*ipd/i to be a probability density and 
ij) to be square integrable on the domain D where our physical system is defined. We 
conclude thus that wave functions representing physical systems must be elements 
of the Hilbert space of square integrable functions defined on the configuration space 
D. We shall denote this space as TL = C 2 (D). The operator H is usually called 
the Hamiltonian operator and it acts on wave functions as a linear differential 
operator. The spectrum will be real, if H is self-adjoint. If it is unbounded and 
if its domain of definition does not coincide with H we say that we have to "face 
domain problems" . The requirement of linearity arises from the superposition rule 
usually associated with the description of interference phenomena. 

To avoid dealing with unbounded operators, often one starts from the "exponen- 
tiated version" of our evolution equation, i.e. a formal solution of the equation [2] 
which we write as 

Tjj(t) = U(t,t )i>(t ). (4) 

If H does not depend on time the expression of the evolution operator U(t, to) may 
be written as: 

U(t,t ) = e-* H( - t - to) (5) 
In this form the evolution operator is bounded and always reversible. It defines a 
one-parameter group of unitary transformations on the Hilbert space defined by the 
wave functions. Unitarity of the evolution means that such evolution is compatible 
with the probabilistic interpretation. Starting with the evolution expressed as a 
one-parameter group of unitary transformations, the operator H, which defines the 
Schrodinger equation, emerges as the infinitesimal generator of the group associated 
with the evolution. 

Thus we have identified the basic ingredients appearing in the Schrodinger pic- 
ture of Quantum Mechanics. It is important to remark the presence of the new 
fundamental constant h in the description of this class of phenomena. It implies 
some fundamentally new aspects of the quantum theory with respect to the clas- 
sical theory. For instance it is known that any measurement process requires an 
exchange of energy (or information) between the object we are measuring and the 
measurement apparatus. The presence of h implies that these exchanges can not be 
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made arbitrarily small and therefore idealized to be negligible. Thus, in the mea- 
surement process, we can not conceive of a sharp separation between the "object" 
and the "apparatus" , so that we can disregard the apparatus altogether. Moreover, 
in the measuring process, there is an inherent ambiguity in the "cut" between what 
we identify as "object" and what we identify as the "apparatus" . 

The problem of measurement in quantum theory is a very deep one and goes 
beyond the scope of these notes. We may simple mention that within the von Neu- 
mann formulation of Quantum Mechanics the "measurement" gives rise to the so 
called "collapse of the wave function" . When we measure some dynamical variable 
(i.e. a physical quantity), the wave function which represents the system is pro- 
jected onto one of the cigcnstatcs of the operator A with a probability that can 
be computed. However we shall not insist on these problem any longer and move 
to the general mathematical structures emerging in the description of quantum 
mechanical systems. 

While we shall deal mainly with finite dimensional Hilbert spaces, we are going 
to consider very briefly the Weyl-Wigner formalism in infinite dimensions and men- 
tion the tomographic formulation of Quantum Mechanics. The reason is that this 
formulation is the most suitable for the study of the quantum-classical transition. 

Before closing this introduction and to better put into perspective the Schrodinger 
and the Hciscnbcrg pictures, we are going to make a few general considerations on 
the minimal mathematical structures required for the description of a physical sys- 
tem. From a minimalist point of view, we identify three main ingredients: 

• a space of states, which we define as S, 

• a space of observables, which we denote as O, 

• a real valued pairing /i:Ox5^1. This pairing produces a real number 
out of a state and an observable, and it is associated with the measurement 
process. To be more precise, /i(A,ip), for A e O and ip e S, will be a 
probability measure on the Borel subsets of M. 

It is not difficult to identify these objects in the framework of classical dynamics 
or of classical statistical mechanics. In Quantum Mechanics we have two prevailing 
pictures: 

• The Schrodinger picture. In this case S is identified with the vectors of a 
Hilbert space H, while the set of dynamical variables (the observables) is 
a derived concept: they are identified with self-adjoint operators acting on 
7i. 

• The Heisenberg picture. Here the situation is complementary: the set of 
dynamical variables (the observables) is the primary concept. Observables 
are assumed to be the real elements of a C*-algebra A. The states are then 
a derived concept and are defined as a proper subset of the space of linear 
functionals on A. 

2.2. The mathematical structure of the Schrodinger picture. In this picture 
the relevant carrier space is a Hilbert space TL. Very often, it is realized as a 
space of square-integrable complex valued functions defined on some spatial domain 
DcR", identified with the configuration space of our system. This configuration 
space turns out to be the configuration space of the physical system that we would 
identify in the quantum-classical transition when considering the so-called classical 
limit defined as H — > 0. 
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States are unitary vectors, which are identified when they differ by multiplication 
by a complex number of modulus one (i.e. an overall phase). This implies an 
equivalence relation on TL: given ip\ and ip2 in TL satisfying = 1 = (t/^IV^), 

we have 

^ = e 4 Vi 9 e [0, 2tt) => Vi ~ ih- (6) 

As a result, physical states are identified with points in the complex projective 
space associated to TL. We shall denote it as VTL or occasionally 1Z(TL) (meaning 
the space of rays of TL). The corresponding equivalence classes will be denoted as 

M - {cj> e TL\4> ~ ^} 

On the other hand, observables are defined as self-adjoint operators acting on 
TL. The pairing between states and observables is defined in terms of the Hermitian 
structure of TL: with any pair ([?/>], ^4) we associate the expectation value function 
eA and defined as 

it is clear that the right hand side depends only on the equivalence class [ip). We 
are using here, as in the rest of the paper, the bra-ket notation of Dirac [13 . 
Evolution is defined on the Hilbert space TL by means of the Scrodinger equation 

ihj t \^)=H\^) (8) 

We can elaborate a little on this equation. Consider now the equation for the 
complex conjugate and transpose of the wave function: 

-ihf t (i>\ = W\H (9) 

where the reality of H has been used. 

Combining both equation, it is immediate to see that, at the level of the Hilbert 
space: 

~{l/>,4>) = -j:({ll>,H4>)-{il>H,4>))\ 

which appears as a continuity equation. 

To be specific we consider the Hilbert space realized in terms of wave functions 
and the Hamiltonian operator as a differential operator associated with a particle 
moving in a potential V(x), we find: 

i^ = -^^ + V(x)^, (10) 
at Zm 

where A is the Laplacian operator defined on D and V(x) is the multiplication 
operator 

V(x)%j) = V(x)i) ip G h 2 (D) 

As we have mentioned, the presence of the configuration space D in the realiza- 
tion of the Hilbert space TL in terms of square integrable functions on D incorporates 
some shadows of the classical limit. As a matter of fact, if we introduce the polar 
form of the wave function as 

i>{x,t) = A(x,t)ei w ^, (11) 
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where A and W are real valued functions and — A 2 — p, we decompose the 
complex linear Schrodinger equation into two nonlinear and coupled differential 
equations. Here the equation for W contains an additional "quantum potential" 

TT ,^ . h 2 AA 

16 dW (gradWH 2 . TT ,^ . 

The other equation takes the form of a continuity equation for the probability 
distribution p: 

^ + divp*=o ff= ^K. (12) 

at m 
The vector field associated with the wave function tp: 

_ dx ihib*Vib -ipVtp* 

v= m = TV~i ' ( 13 ) 

at m ip*ip 

is the extra ingredient of Bohmian mechanics ([5]). 

Let us indulge now on the integration of the Schrodinger equation. When the 
Hamiltonian operator does not depend explicitly on time, i.e. it is time-translation 
invariant, we may use the separation of variables and consider 

ip(x,t) = <j)(x)e-^ Et . 

We thus get the Helmholtz equation associated with the Schrodinger equation: 

h 2 -\ 
-— A + VU = ^ 
2m / 

One solves for this equation in terms of an "eigenvalue problem" with preassigned 
initial conditions (boundary conditions) and finds, for suitable potentials, a set of 
fundamental solutions which provide also a basis for the Hilbert space, say 



with the corresponding eigenvalues 

Ei , Ei , • • • 

A general solution for our quantum problem will thus have the form: 

^x,t) = J2^M^ Eit , (14) 

3 

with ij}(x, 0) = ■ Cj(f>j(x) allowing to determine {a,} up to an overall phase. For 
simplicity we have considered a situation where eigenvalues are not degenerate. 
The merit of the particular basis we have considered relies on the fact that on 
these vectors H acts as a multiplication operator, i.e. as a diagonal matrix. Let 
us elaborate a little more on this aspect in more general terms. To this aim it is 
convenient to go back to the abstract Dirac notation in terms of bras and kets. 
Thus starting with Schrodinger equation 

> = >■ 

Let us consider now an orthonormal basis for the Hilbert space {le;)}^^... and the 
corresponding basis for the dual space {(ej|}i=i,..., with the property 

(ej\ei) = 5 jk 
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In this basis the Hcrmitian structure is represented by the identity matrix. The 
orthogonality also allows to write the decomposition of the identity operator 



With the help of these bases we can write our initial Schrodinger equation in the 
matrix form: 

3 jk 

Now, by using the fact that {|ej)} are a basis, we get: 

k 

Denoting h*- = X)fc( e ilj|l e fc) we can write the Schrodinger equation as 

^ = 

k 

This equation would allow us to use all the results from matrix algebra if the 
matrices would be finite dimensional. As a matter of fact, this would be the case if 
for any chosen j the number of indices k coming into h° k were finite. In geometrical 
terms this would mean that we can decompose the full Hilbert space into the 
direct sum of invariant finite dimensional Hilbert spaces. In this way the infinite 
dimensional matrix associated with the operator would decompose into a block 
diagonal form, each block being a finite dimensional matrix. Thus the eigenvalue 
problem associated with the Hclmholtz equation is nothing but a way to find such 
an invariant decomposition of the Hilbert space where, in addition, each finite 
dimensional block is a multiple of the identity, the multiplication factor being the 
eigenvalue of the corresponding eigenspace. This procedure gives rise to the well 
known spectral decomposition of the Hamiltonian operator. 

Very often, an invariant decomposition of the full Hilbert space will be associ- 
ated with the decomposition into irreducible representations of a compact group of 
symmetries for the Hamiltonian operator. 

In conclusion, the most convenient basis {|ei)} depends on the particular problem 
we are dealing with, and there are no general prescriptions. 

Of course, due to the role of the observer, or apparatus, or laboratory, a conve- 
nient basis from the point of view of interpretation, is provided by the eigenstates 

of the particle's position vector operator X (whose components are the three com- 
muting coordinate operators xi,X2,x$, defined as Xj\x) = Xj\x)). 

This continuous basis satisfies some generalized orthonormality conditions, namely 



\x)d 3 x{x\ = I (x\x") = S(x-x'). (15) 
Decomposing \<fi) in the \x) basis yields 

\x)d 3 x(x\<p). 



/' 



Usually we write (x\(j>) — <p(x) and this basis provides us with a clear interpreta- 
tion of the wave function, i.e. the component of the abstract vector \(p} in the basis 



BASICS OF QUANTUM MECHANICS, GEOMETRIZATION AND SOME APPLICATIONS ... 9 



of the position operators. Clearly the continuous index x replaces previous discrete 
indices k 7 j and integration replaces the summation. 

The position representation is quite suitable for measurements of the particle's 
position. In other instances, as for instance the case of massless particles, the 
momentum representation is more convenient. In this case we use the momentum 

operator basis \p) associated with the momentum operator P: 

P\P) =P\P)- 

Again a decomposition of the identity can be written with this basis: 

J \p)d 3 p(p\ = I (pV) = S(P- PO- (16) 
Each of these states can be decomposed in the position basis: 

\P) = J \x)d 3 x(m = j^J^ J d 3 xeinx) 

Thus the components of the momentum eigenstates, in the position representa- 
tion, are given by plane-waves. Therefore they are not square-integrable functions 
and can not be considered as elements of L 2 (D). In this formalism, to make sense 
of the eigenstates of the position or of the momentum operators, one usually deals 
with wave packets. These are usually written as: 

= (2^^ / d 3 pcp{p)^ S - Et \ (17) 

and E is written in terms of the momenta. The particular expression of E in terms 
of the momenta is usually called the dispersion relation and the integral appear- 
ing in the momentum variables should be considered as taken on the submanifold 
defined by the dispersion relation. 

The Fourier transform which appears in going from one set (of maximally com- 
muting) observables to a conjugate set (of maximally commuting) variables plays a 
very important role because it is a unitary transformation and therefore compatible 
with the probabilistic interpretation. 

Moreover, from the mathematical point of view, we notice that the Fourier trans- 
form represents the "harmonic analysis" on the Abelian vector group. Thus it is 
associated with translations in space-time or translations in the energy-momentum 
space [30]. 

It is interesting to write now the evolution operator we mentioned above Q, 
when we consider the position representation. We find: 

(xm)) = J d 3 x'{x\U{t,t )\x'){x'\i>{t )), 

i.e. 

t/;(x*,t) = J d 3 x'G(x,x";t,t )il){x",t ). 

Here we have introduced the standard notation for the Green function, the "prop- 
agator" which makes our initial state to evolve. It simply represents the matrix 
form of the evolution operator in the continuous basis \x) (in particular the matrix 
element (x\U(t,to)\x')). It is simple to see that at the initial time 

G(x,x';0) = S(x-x'). 
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Therefore these matrix elements become functions on Q x Q (Q being the configu- 
ration space). By using J d 3 x\x}{x\ or j d 3 p\p)(p\, we can find analogous matrix el- 
ements written as functions of the cotangent bundle T*Q, i.e. G(p, x; t) or G(x,p; t) 
respectively. 

The probabilistic interpretation we have considered in the introduction is now 
formalized by saying that when the state of our system is described by the wave 
function ijj(x,t), the probability of finding the particle in the region D C M. 3 at time 
t, is given by 

T(D;t)= [ \iP(x,t)\ 2 d 3 x 

The operator 

V D = ( \x)d 3 x(x\ (18) 
Jd 

is a projection operator and may also be described by means of the characteristic 
function 

"1 if x e D 
ifx(£D 



Xd{x) 
Thus we can also write: 



V D = J XD(x)\x)d 3 x(x\. 



Similar formulae can also be written for the momentum representation or any 
other realization of the Hilbert space in terms of square integrablc functions. 

Up to here, we have been using the position operator representation or the 
momentum-operator representation and mentioned the possibility of using either 
one. One may wonder if the general aspects could be dealt with without making 
an a-priori choice of the representation. 

As a matter of fact, it is indeed possible and this is what is done in the Weyl- 
Wigncr formalism along with the further elaboration given by the tomographic 
formalism. 

The reduction to a specific representation is due to the requirement of irreducibil- 
ity of the representation of the canonical commutation relations. In the following 
sections we shall briefly consider the Weyl-Wigncr formalism, the canonical com- 
mutation relations and the emergence of these irreducible representations. 



3. Weyl systems 

3.1. Definition and main properties. A Weyl system is defined on a symplectic 
vector space (V,lj) in the following way. It is a map 

W: (V,w) -» U(H), v^W(v), (19) 

from the symplectic vector space to unitary operators U(H) on some Hilbert space 
Ti. We require that this map defines a strongly continuous one-parameter group of 
transformations W(tv) and that it satisfies the property 

W{v!)W(v 2 ) = e^" 1 '" 2 ^^! + v 2 ) y Vl v 2 e V 

Strong continuity, by means of Stone's theorem, implies the existence of an 
(essentially) selfadjoint operator R(v) such that 

W(tv) = e itR ^ 
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A well-known theorem by von Neumann establishes that for any finite-dimensional 
symplectic vector space a Weyl system can always be obtained in terms of operators 
acting on square-integrable functions defined on any Lagrangian subspace L <Z V , 
with the translational invariant Lebesgue measure We recall that a Lagrangian 
subspace of a symplectic vector space is a maximal vector subspace L such that 
w{y\, V2) = for any pair of vectors t>i,i>2 £ L. Maximality implies that any other 
vector u £ V such that u>(v, u) — for any d£L must belong to L. 

It is immediately seen that when vectors Ui,i>2 are such that u){v\,V2) — 0, the 
operators W(v\) and W(y2) commute. Therefore, a Lagrangian subspace arises as 
the joint spectrum of a maximal commuting set of unitary operators. 

By using vectors \y), eigenvectors of the operators associated with vectors of L, 
we can write a decomposition of the identity 

1 = J \v)Hl{v\, 

and also 

\ip) = / \y)nL(y\ip) = / ^(y)\y)nL- 

J L J L 

Within this setting, a Weyl system can be constructed in the following way. We 
consider the Lagrangian subspace L c V and its cotangent bundle T*L~Lxi*. 
which is endowed with the canonical symplectic form u>q. Now we construct a 
symplcctomorphism 

(v,w) — (r*i |Wo ). 

Thus we can write the elements of V as pairs v — (x, a) where x £ L and a £ 
L* . With this decomposition, Weyl operators will be defined by their action on 
L 2 (L,fi L ) as 

[W{x,0)i>}{y)=i>{x + y) x,y £ L t/j £ L 2 (L, fi L ) 

[l¥(0,a)V](2/) = e M ^V(y) y£L,a£L* ^£L 2 (L,^ L ). 
We know from Stone's theorem that there exist infinitesimal generators for these 
transformations. We may introduce the operators 

J« =#([,;()) Qj = i?(0; 

where we use lj to represent the base vectors for L and L* , i.e. lj = (0, 0, • • • , 1, 0, • ■ ■ , 
1 being in the j-th position. 

With these operators, the Weyl map has the form: 

W(x,a)=eM w * pi+a '4'\ (20) 
Equivalently we could write instead: 

Wi(x,a) = e^^es*^, 

or 

W 2 (x,a) = ei x i p3 ei a3( i>. 

These various forms differ by multiplication by a complex number of modulus 
one and give rise to different "orderings" (see [1]). 

If we select a complex structure on V, defined by a tensor J : V — > V satisfying 
J 2 = —I, we can define two operators 
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a(v)=^(R(v)+iR(Jv)) a+(v)==^(R(v)-iR(Jv)), (21) 

which will be called annihilation and creation operators respectively. Also the 
natural symplectic structure u — da? A dxj can be written in the complex coordi- 
nates 

zj = (Xj + ia 1 ) z* = -j= ( Xj - ia> ) 

as 

u> — idz* A dzj . 

Trivially it follows that (z\, ■ ■ ■ , z n , ■ ■ ■ ) are coordinates of a Lagrangian subspace, 
as well as (z* , ■ ■ ■ , z* , • • • ) . 

In these coordinates, we have the Weyl map given by 

z i ► D{z) = e za+ - z * a . 

In this form, Weyl operators are called displacement operators and are used to 
construct "coherent states" . 

By means of a symplectic linear transformation, T : V — > V such that ui{Tv\, TV2) 
oj(vi, V2), we may define an automorphism on the set of unitary operators by setting: 

W(Tv) = v T (W(v)) = U(T)W{v)U\T), 

because automorphisms of the unitary group are inner. 

Thus with a one-parameter group of linear homogeneous symplectic transforma- 
tions, we define a one-parameter group of unitary transformations. As symplectic 
linear transformations have an infinitesimal generator associated with a quadratic 
function on V, the corresponding infinitesimal generator of the unitary transforma- 
tion (Stone's theorem) shall be the operator associated with the quadratic function 
on V. In conclusion we obtain: 

Lemma 1. We can associate self-adjoint operators with quadratic functions on V . 

By using the set of complex coordinates introduced above, and the association 

z 1 ► W(z) = e za+ ~ z " a , 

we can construct the automorphism 

A 1 ► A(z) = W\z)AW{z). 

By using the vacuum state |0) as a fiducial vector, defined for instance as a|0) = 0, 
we may set a one-to-one correspondence between Hermitian operators and functions 
on V by setting: 

f A (z) = (0\A(z)\0) = (z\A\z). 

This function is known as the Berezin symbol of the operator A and allows to 
define an invertible mapping between Hilbert-Schmidt operators and square inte- 
grable functions on V. 

We should also remark that the map 

z 1 ^ W(z)\0) = \z) 

is a way to immerse C™ (the phase space) , as a submanifold (or a subset) of the 
Hilbcrt space H (say, the Fock space). But it is not a linear map onto the image. 
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If we consider the expectation value functions 

(m\i>) 

we may pull-back them to C" by considering: 

(z\A\z) 

e A (z) - 



(z\z) ■ 

Another important point regards the algebraic structures on the sets. As the 
mapping (24) is one-to-one, it is possible to associate a function to the product 
of two operators, and thus define a binary operation on functions associated with 
operators by setting 

A^ f A (^),B^ f B (^) AB^f AB (ij;) = (iP\AB\iP) = {f A *f B )(iP). (22) 

When we pull-back this product to the phase space (i.e. define (f A *f B )(z)), we 
get a star-product very close to the usual Moyal product ([30]). As a matter of 
fact 

J lm n n = IJaJb}, 

ft— *0 ti 

and thus we can recover the Poisson bracket (associated to the symplectic form u>) 
on phase space. 

Let us summarize then the results presented so far: 

• Weyl maps are applications from a symplectic vector space ( V, u>) to the set 
of unitary operators of some Hilbert space U (Ti.) satisfying two properties: 

— the map is strongly continuous, 

- W(v!)W(v 2 ) = W(v 2 )W(vi)eT; u{ - v i- v ^ for any v u v 2 G V. 

• Stone's theorem allows us to define the infinitesimal generator of the unitary 
transformation corresponding to the vector v. We denote this infinitesimal 
generator as iR(v): 

W(v) = e iR ^ Vv G V. 
The condition above implies that 

[R(v{),R(v2)] = ihw(vi,v 2 ) V«i,«2 G V 

• If we realize the symplectic vector space V as the cotangent bundle (with 
its canonical symplectic form) of some Lagrangian submanifold L a V, i.e. 
V ~ T* L ; we can write the expression of the Weyl map as: 

W(x, a) = e^ xP+a ® (x, a) G T*L; 

where the operators P and Q via Stone's theorem and have an associated 
action on the elements of L 2 (L, i^i) as 

(Qi>)(y) = yi>(y) (P^)(y) = -it&- i> eI?(L^ L ). 

oy 

• We can also introduce the set of complex coordinates (z, z*) and write the 
expression of the Weyl map as 

W(z) = e za+ - z * a ; 



for a and a + the creation and annihilation operators (21). 
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3.2. Weyl maps and Poisson tensors. Let us elaborate a little further and 
consider a similar construction, but having as initial set the Poisson vector space 
(V*,A), i.e. the dual space V* — Lin(V, R) and the Poisson tensor A defined from 
the symplectic form oj as: 

A(a Vl ,a V2 ) = w(ui,u 2 ), 

for a Vi the element of V* defined as 

a Vi {u) = uj{vi, u) Vu G V. 

With these elements, we define the Weyl map as an application 

W: (V*,A)->H, 

being strongly continuous, and satisfying the condition 

W{ ai )W{a 2 ) = W(a 2 )W(a 1 )e lHai ' a2) . 

If we use this as our starting point, we may define Weyl systems also for de- 
generate Poisson structures. This turns out to be very useful when we deal with 
constrained systems in the sense of Dirac and Bergmann |15j . 

3.3. Weyl systems and linear transformations. Consider again the symplectic 
vector space (V, ui) and a linear transformation T : V — > V. We can define a new 
tensor: 

ujt(vi,v 2 ) = lu(Tvi,Tv 2 ). 

If T is invertible, lot is again a symplectic form. We can thus associate a Weyl 
system to the new symplectic form lot, in the form: 

(V,w) ^ >U{H) 



(V> T ) — -+U(H) 

As the diagram is commutative, vt(W(v)) — Wt(Tv). 

At the infinitesimal level, we have that Rt(v) = R(Tv) and therefore Rt(v) is 
not Hermitian with respect to the Hermitian structure on H which makes W(v) 
unitary. Hence, this is a situation where Rt(v) may have a real spectrum without 
being Hermitian. This situation has been considered in the literature as the one 
using pseudo-Hermitian operators ([2| I4T1 129] ). 

When T is a symplectic transformation, and hence = to, the transformation 
vt is an inner automorphism: 

v T (W{v)) = U T W{v)U T . 

If we consider now a one-parameter family of linear symplectic transformations 
T(A), where A G M; we know that the infinitesimal generator of this family is a 
Hamiltonian vector field. As the transformation is linear, we also know that the 
associated Hamiltonian function is quadratic. But we can also consider the infini- 
tesimal generator of the associated unitary transformation Ut(\)- This infinitesimal 
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generator will be a self-adjoint operator on H, which is thus associated to the qua- 
dratic Hamiltonian function on V which generated the family of linear symplcctic 
transformations . 

As an example, if we consider the transformation 



T = 



-1 

1 



the associated unitary operator is given by the Fourier transform. 

Let us consider now the Heisenberg-Weyl Lie algebra (V © R, u) defined as 

[vi,v 2 ] = u>(vi,v 2 ); 

and its holomorph, i.e. the extension of the algebra with its derivation algebra. Let 
us recall that the derivation algebra of a Lie algebra is the set of linear transforma- 
tions A : V — > V which satisfy: 

A([v 1 ,v 2 }) = [A(v 1 ),v 2 ] + [v 1 ,A(v 2 )}. 

In the present situation, this corresponds to the algebra of infinitesimal canonical 
transformations of V. If we consider then this extension of the Hciscnbcrg-Weyl 
algebra, the Weyl construction shows that the entire associated Lie group can be 
represented by unitary operators on the Hilbcrt space on which the Weyl operators 
act. 

3.4. Associating operators to more general functions. We have seen that 
our previous construction is able to provide a unitary operator to represent any 
function of V which is at most quadratic. What can be said for more general 
functions? 

The trick elaborated by Weyl was the following. Consider the realization V ~ 
T*L which assigns a pair (x,a) to any element v £ V. Now, we can consider the 
double Fourier transform of a function / : V — > R as 



f(q,p) = J J dxdae^+^f(x,a). 

Now define the operator / associated to the function / as 

f{Q,P) = J dxdae^^ +xp \f(x,a). (23) 

This map can be inverted, i.e. we can find a map from operators acting on Ti 
onto functions defined on V, such that it is an involutive map of index two. The 
inverse, namely the function Ja associated to an operator A on the Hilbert space, 
can be written as: 

f A (v) = TrAW(v) (24) 
We might also use the symplectic Fourier transform 

f(q,p) = J dxdae-'^^^^TviAw^x^)). 

Definition 1. If we apply this to the density operator associated to a pure state 
\ip), — \ip)(ip\, we obtain: 

This is called the Wigner function associated to the state \ip). 
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We can conclude thus with the following statement: 

Lemma 2. The Wigner-Weyl map defines a bijection between Hilbert- Schmidt op- 
erators and square-integrable functions on phase space. 

If we also denote by W the map / i— ► / = W(f) (i.e. applied to functions on 
V), we can use the bijection to transport the associative product of the C*-Weyl 
algebra generated by unitary operators on 7i onto the set of functions on V. This 
defines the ^—product we saw above: 

(f*g)(x,a) = W- x {W{f).W{g)){x,a) = Tr (W(f)W(g)W*(x, a)) . 

This is the integral form of the Moyal product, which is usually written in terms 
of bidifferential operators 

(/*$)(«,*) = /(^^"^Sfep), (26) 

or equivalcntly 

(f*g){q,p) = /te + ^jL'P-^J^sfozO- 

It should be clear now that all descriptions in terms of operators have a coun- 
terpart in terms of complex valued functions on the phase space with the Moyal 
product. For instance, we can write the eigenvalue problem for the Hamiltonian of 
our quantum system as: 

H-k p E = Ep E , 

where H is the function corresponding to the Hamiltonian operator and pe is the 
function representing the cigenstate with energy E. Similarly we can write evolution 
equations and their exponentiation as a Taylor expansion in terms of ^—products. 

3.5. Tomograms. The Wigner function associated to a state ip was introduced 
in 1932 as a tool for the study of quantum corrections to classical equilibrium 
distributions [32]. However, unlike the distribution function of classical statistical 
mechanics which is a probability distribution function on phase space, the Wigner 
function is only a quasi-distribution, i.e. it is not positive definite. 

The tomogram description of quantum mechanics uses a formulation in which the 
quantum state is described by conventional non-negative probability distributions. 
Thus, within this formulation both the classical and the quantum descriptions of a 
physical system use true probability distributions, and this property makes easier to 
compare predictions of the two descriptions. Moreover, it makes also easier to deal 
with the quantum-classical transition because the describing functions are defined 
on the same carrier space. 

We shall not enter in this fascinating subject, but instead refer to some previous 
work to find details on this formulation in the same spirit of these notes [231 H3] 

4. The geometrical description of Quantum Mechanics 

4.1. The Hermitian structure. To consider Quantum Mechanics from a differ- 
ential geometric point of view, it is convenient to consider the Hilbert space 7i as 
a real differential manifold. In order to do this, we have to convert the Hermitian 
inner product (•, •) which defines TL into an Hermitian tensor. This is quite similar 
to what we do in the transition from special relativity to general relativity. In that 
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case, the Minkowski product r] ilv x ll x v is replaced by the tensor -q^dx^ ® dx v and 
thus space-time becomes a manifold whose tangent space becomes endowed with 
a Minkowskian product. In the present case, the Hilbert space H is replaced by a 
Hilbertian manifold, whose tangent space is a Hilbert space with the inner product 
(;■)■ 

To be definite, let us consider an orthonormal basis {|ej)}j=i, -- for H and define 
the coordinate functions: 

(ej\ip) = z J (ip) = q 3 (ip) +ipj(ip), 

where q 3 and pj are real functions. With this in mind, we replace the inner product 
(ip\ip) by (dip ® dip), where dip is a Hilbert-space valued one-form, i.e. a map from 
TJi to Ji linear along the fibers: 

^| ej )(e,|#)=^^| ej ) 

j j 

Then the expression of the product becomes: 

(dip <g> dtp) = dz 3 (ej\e k }dz k = dz 3 <g> dz k (e j \e k ) 
Now using the real functions q 3 ,Pj introduced above we get: 

(dip ® dip) = (dq 3 - idpj) ® (dg fe + zdpfc) (ej |efe) = 

(dq 3 ® + dpj <g dp/c + i(cV <g> dpk — dpj (g) dq h )) (e \ek) 

Because of the orthonormality of the base we find: 

(dip <g> = (dq 3 <g> <ig fe + dp 3 <g> dp fe + i(dtf* <g> dp k - dpj <g> dq k )) S jk 

We can also consider the scalar product on the dual space H* and thus introduce 
two tensors: 

d d d d 
dq k dq k dp k dp k 



~*ih f)n, /-)<n& ' 



and 



Clearly, 



d d d d 
dpk dq k dq k dp k 



Op k Oq k dq k dp k 
defines a Poisson braket on the space of functions while 

G df, dg) ^ fc <g> ^ + ^ fc <g> ^ 

defines a commutative bracket. 

And we can define the Hermitian bracket as a combination of the above products: 

(df,dg) = G(df,dg) + iA(df,dg). 
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J = E U'aor-^sT]. (29) 



4.2. Other geometrical objects. We may use the set of real coordinate functions 
introduced in the previous section to write explicitly the complex structure of Ti. 
We know that the Hilbert space is complex, therefore if we want to describe it 
as a real differentiable manifold, we need to identify the corresponding complex 
structure. It corresponds to a (l,l)-tensor J, satisfying the property J 2 = —I, 
which written in the coordinates above has the form 

d , d 

dp k <gs — - 

The other ingredient of the Hilbert space structure of Ti is its linear structure. 
This also can be encoded in a tensor, a vector field in this case, which is the dilation 
vector field A : Ti — ► Ti x Ti ~ TTi. It is defined as A(ip) — (ip, ip) and hence takes 
the coordinate expression: 

The complex and the linear structures can be thus combined as 

TV AN ® k ® 

which is the infinitesimal generator of the multiplication by a phase. 

Let us go back now to the problem of the star product at the level of the functions. 
Consider an operator A on Ti and the evaluation functions 

f A (il>) = (1>\A\il>), 

and the expectation value function 



ex WO 



w#> 



The product which translates the associative algebra structure of the set of 
operators reads 



and 



1 i 
fA*fB = ^G(df A , df B ) + -K{df A ,df B ) = f AB . 



1 i 

e A *e B = e A e B + -G P {de Al de B ) + -K P (de A , de B ) = e AB , 



where Gp = (A|A)G and Ap = (A|A) (we shall go back to these objects in the 
next sections). 

Definition 2. A Kahlerian function is a function f such that the Hamiltonian 
vector field Tf, A(df) = Tt, is also a Killing vector field, i.e. 

Lv f G = 

On Kahlerian functions, G defines a Jordan algebra structure, while A defines 
a Lie structure. Both structures are compatible (they define a structure called 
Lie- Jordan algebra [HI [22]) an d together define the C*-algebra corresponding to 
bounded operators. 

The operator norm can also be obtained from the function as \\A\\ = swpe A {ijj) 
for any positive self-adjoint operator A and any normalized vector tp. 
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If A is a general, possibly non self-adjoint operator, we have that \\A\\ 2 = ||j4*v1|| 
and this is translated as 

\\A\\ 2 = sup^(f A *f A )W 

Remark 1. The functional representation of a C* -algebra requires a uniform 
Kahler bundle |10l [9, 8 p : V — > B. V will be the set of pure states of A, while B 
is the spectrum of A (i.e. the set of unitary equivalent classes of irreducible Hilber- 
tian representations of A), p is the natural projection, associating the equivalence 
class of irreducible representations containing the GNS representation induced by 
the pure state co: 

(J^(A) = e n ( A )(ip) = Tray7r(A) \/A e A 

For further details on the geometrical interpretation of the GNS construction see 
Chruchinski-Marmo [7] . 

Let us consider now the relation with the notion of Gclfand transform. The 
Gelfand transform of an operator A is defined as 

e A :T -> C V 3 lo^ e A (uj) = uj(A) 

The Gelfand transform A <— > e A is a linear involutive preserving injection of A into 
T{V). For A, B € A, e A B — e A *eB and 

\\ A \\ 2 = sup U £ V (e A *e A )(u)- 

The range of the transform is the set 

K-uiV) = {/€ JC\f, /*/>/* I are uniformly continuous on V} 

4.3. Projective Hilbert spaces as Kahler manifolds. We already saw above 
that from a physical point of view, the probabilistic interpretation requires that the 
set of states of a system is not a vector space but a complex projective one. Our 
aim now is to present the geometrical structures arising in this case. 

First let us study the expression of the Hermitian tensor in VTi. It is simple 
to obtain a Hermitian tensor field on Ti which vanishes on vertical vector fields 
(generated by A and J (A)) and provides a projectable function when evaluated on 
projectable vector fields: 

{dip® dip) (ip\dtp) ® (dtp\tp) 

From a geometrical point of view, the complex projective space is obtained as the 
quotient manifold obtained from the Hilbert space TL with respect to the foliation 
generated by two vector fields: 

• the dilation vector field A = q k + pk g^- , and 

• the phase-change vector field obtained as the vector field resulting from 
the action of the complex structure J on A, i.e. T = J(A) = 9 fc g^ — 
Pk -g^k ■ This vector field is the Hamiltonian vector field (with respect to the 
canonical Poisson structure A of the function f(q,p) — (ip\ip), associated 
with the unity operator I.. 

From this observation it is simple now to study, from a geometrical point of 
view, the projection of the contravariant tensors introduced on TL to the complex 
projective space. We need to consider contravariant tensors and functions which 
are invariant under the action of T and A . 
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From this perspective, it is simple to see that the functions that we have 
associated with operators A are trivially projectable. On the other hand, it is 
simple to see (from the coordinate expressions, for instance) that the tensors G 
and A are homogeneous of degree -2, and hence non-projectable. The homogeneity 
consideration is useful, though, to modify them by a conformal factor and make 
them projectable. Thus the objects 

G p = (A|A)G A p = (A|A)A, (31) 

define projectable tensors. 

Particularly the issue on the Poisson tensor has some deep implications from a 
geometrical point of view. 

We may use the two tensors above to realize the C*-algebra of operators on Ti as 
a C*-algebra of functions on the complex projective space. With the help of these 
functions we can formulate the eigenvalue problem on the complex projective space 
considered as a real manifold, and write the equations of motion corresponding to 
the Heisenberg formulation of the evolution equations. Thanks to the Wigner's 
theorem, it is possible to solve the equations of the motion on the Hilbert space, 
by exponentiation for instance, and then to project the result onto the ray space 
to obtain the desired solutions. More likely it is because of this property that 
physicists have not bothered dealing with Quantum Mechanics on the ray space 
considered as a real manifold. 

4.4. The momentum mapping. To end this quick description of the basics of 
the geometrical description of Quantum Mechanics, we will now consider another 
important aspect, here worked out for a specific example. 

Given the Hilbert space Tt = C 3 (or equivalcntly the corresponding projective 
space VH), we know that there is an action of the unitary group U(3) on them. 
With respect to the canonical symplectic structure of H this action is strongly 
symplectic (|26j). and therefore admits a momentum map 

H:H = H- {0} ^u*(3) 

or equivalcntly 

These mappings commute with respect to the natural projection ir : H — ► VH, 

as 



Ho 




vn 



Given a point 6 Ti, its image by fi is an element of the dual of the algebra of 
operators. From that point of view, it is easy to construct a pairing between states 
and operators defined as 



Hty)(A) = (i>\A\i/>) = p+{A) = Tr(A Pi> ) = Tr(A^)(^|) (32) 
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Thus wc define the momentum mappings as 

/i(V0=/ty = MM pCP) = Pm = l ^M, (33) 

where IV'KV'I i s the mapping defined on TL as 

\ip)(ip\:H^H H ^4>^ (tp\(f>) \if>). 

Clearly the mapping p is one-to-one. This implies that the geometric structures 
that we know on VTi must have an analogue on u*(3). Indeed, it is simple to see 
that the geometric structures correspond to the canonical structures of the coadjoint 
orbit of J7(3) on u*(3) which contains p([ip)). 

One important fact is the isomorphism u(3) — > u*(3) given by the Killing-Cartan 
metric. Indeed, there is a one-to-one correspondence between the elements of the Lie 
algebra and the elements of its dual and this correspondence intertwines the adjoint 
and the coadjoint actions. We shall denote as a a the element of u*(3) corresponding 
to the observable iA for A £ u(3). Besides, as the set of physical observables in this 
case corresponds to the set of linear functions on u*(3)(u(7i) = Lin(u*(7i), M.)), we 
can use the pullback p* : JF(u*(3)) — ► J-{Ti) to relate the set of Hermitian operators 
and the set of functions defined on the Hilbcrt space. The geometrical structures 
are also related, defining suitable morphisms: 

Proposition 1. i) The momentum map is equivariant with respect to the ac- 

tion ofU(H) on Ho and the coadjoint action of U(H.) on u*(H). In par- 
ticular, this says that the Schrodinger equation of motion on H is p-related 
with the Heisenberg equation of motion on u(7i) (the space u(H) is iden- 
tified with the dual by means of the scalar product defined by the trace). 
Moreover, 

ii) n*{A) = f A , p*{A) = e A . 

iii) p*({A,B}) = and p*({A,B}) = {e A ,e B }, where 

iv) p* (R(dA, dB)) = G(fi* (dA), fi* (dB)) and for the other mapping 

p,*{R{dA,dB)) = G P (il*(dA),il*(dB)) + e A e B , 

where R is the Jordan tensor defined on u*(3) from the canonical one on 
u(3) as 

K(0(dA, dB) = f([A, B]+) = '-T^AB + BA) , (34) 

for £ S u*(3) ; and A and B are arbitrary elements o/u*(3). 

Proof. Direct computation. 

Instead of insisting with additional general aspects of the geometrical formulation 
of Quantum Mechanics (further details may be found in [5J |5] , we shall consider 
now some specific examples which will illustrate the general picture. 

4.5. States: Density states. We have seen how the momentum map fx allows us 
to embed the complex projective space PTt on the dual of the Lie algebra u(TC). 
The resulting elements represent the set of pure states of the quantum system. 
But in many physical situations we have more general states, i.e. density states 
which are convex combinations of pure states. They are represented by a family 
p = {pi, ■ ■ ■ , pt\, each element satisfying 

Pk = Pk, Pk = Pk, Trp fe = 1, 
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along with a probability vector, namely p — (pi,P2) ■ ■ ■ ,Pk) with ^2jPj — 1 and 
Pj > Vj. Out of these we construct a density state p = ■ PjPj. The evaluation 
of this state on some observable A is given by 

p{A)=Y jPj Trp j A = TrpA. (35) 

3 

We shall call density states all convex combinations of pure states, and we 
denote them by V{H) [SHE]. 

As any of the elements in p can be embedded into u*(7i), it makes perfect sense 
to consider p also as an element in the dual of the unitary algebra. And we hence 
consider the geometric structure we defined on u* (TL) as the Poisson or the Jordan 
brackets 

{fA,fB}(p) = ^2Pkf[A,B]_(lpk) = ^2Pk{fA,fB}(^k) 
k k 

(f A J B )(p) = J2Pk(fAjB)(^k) (36) 

k 

where f A (p) = EjfePfeAW' 

As for the geometric structures on T>(H) we shall consider it as a real manifold 
with boundary embedded into the real vector space u*(Tt). On this space the two 



structures above (36 1, define a Poisson and a "Riemannian" structure. The Poisson 
structure is degenerate. However it is also possible to define a generalized complex 
structure satisfying 

J 3 = -J (37) 
The boundary is a stratified manifold, corresponding to the union of symplectic 
orbits of U(H.) of different dimensions, passing through density matrices of not 
maximal rank. For further information see [201 [ST1 155] , 

5. Example I: two level quantum systems 

5.1. The Hilbert space. For a two levels system we shall consider an orthonormal 
basis on C 2 , say {|ei), |e2)}. We introduce thus a set of coordinates 

(ej\i/>) = z J (V0 = q 3 W + iPjW 3 = h 2 - 
In the following we shall use z 3 or q° , pj omitting the dependence in the state ip 

as it is usually done in differential geometry. 

The set of physical states is not equal to C 2 , since we have to consider the 

equivalence relation given by the multiplication by a complex number i.e. 

i)x ~ ^ V-2 = A^i A S C = C - {0}. 

And besides, the norm of the state must be equal to one. These two properties can 
be encoded in the following diagram: 

C 2 -S 2 



s 3 

where S 2 and S 3 stand for the two and three dimensional spheres, and the projection 
th is the Hopf fibration. The projection ir is associating each vector with the 
one-dimensional complex vector space to which it belongs. Thus we see how this 
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projection factorizes through a projection onto S 3 and a further projection given 
by the Hopf fibration (a [/(l)-fibration). 

The Hermitian inner product on C 2 can be written in the coordinates Z\,Z2 as 

(ip\ip) = z j z k (e k \e j ) = ZjZ . 

Equivalcntly we can write it in real coordinates q,p and obtain: 

=pi+p 2 2 + ql + <fi 

The associated tensor field reads 

(dtp <g> dtp) = dz~j <g> dz (efe|ej), 

or in real coordinates 

(dtp ® eh/;) = (rf^j — idp?) <g> (rfg* 1 + idpfe)(e/;|ej). 
This expression can be decomposed into its real and its imaginary parts and obtain 
dq k ® dq k + dp k <g> dp/j tig, ® tip' — dp? ® dqj = dq rj A dp? . 

Hence we recognize a Riemannian structure, on the left (the real part) and a 
symplcctic structure (the imaginary one, on the right). 

We can also obtain these tensors in contravariant form if we take as starting 
point the Hilbert space (C 2 )*. If we repeat the steps above, we obtain the two 
contravariant tensors: 

_ d d d d d d 

dq k dq k dp k dp k oq k ap k 

Other tensors encode the complex vector space structure of H = C 2 : 

• the dilation vector field A = q x ^ + p 1 ^ + q 2 ^ + p 2 ^ , 

• and the complex structure tensor J = dp\ ® ^ — dq 1 ® -J^ + dp 2 ® ^ — 
dq 2 ® gfi- 

By combining both tensors, we can define the infinitesimal generator of the 
multiplication by a phase: 

t.wan d , d d 9 d 

r = J(A) =ft- 9 7TT+P 2 ^ 

aqi op L oq 2 op z 

Thus we see how A is responsible for the quotienting from C 2 onto S 3 , while T 
is responsible for the Hopf fibration S 3 — > S* 2 . 

By using the Hermitian operators {<7o, cti, 02, 03} to construct functions (i/)|A|i/>), 
we obtain the real quadratic functions 

<7i +p? + <?2 +pi qm+PiP2, qiP2-piq2, ql + p\-{ql+pl)- 

It is not difficult now to compute the Poisson brackets of these quadratic func- 
tions to find that they arc the Hamiltonian for the infinitesimal generators of the 
u(2) algebra. We may also compute explicitly the Jordan brackets, as for instance 

{q\+p\ + q\ +pl,qiq2+PiP2}+ = 4(<7i<7 2 + PiP2)- 

Similar results are obtained with the other brackets. The result we want to point 
out is 

Lemma 3. The function q 2 +p\ + q\ +p\ acts, with respect to the Jordan bracket, 
as the identity operator except for a normalization factor. 
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We also find 

Um +PiP2,qiq2 + Pip 2 }+ = 4(<?i + p\ + ql+p%). 
And analogously for the other quadratic functions. We have 

Lemma 4. The product of all functions (in the family above) with themselves 
produce a multiple of the quadratic isotropic function. 

We also can obtain easily 

{q\q2 +piP2,qiP2 -pm}+ = = {q x q 2 + piP2,ql +p\ - {q\ +p\)}+ 

The additional relevant property is that the Hamiltonian vector field are Killing 
vectors. In terms of brackets this amounts to: 



This condition, plus the compatibility between the associator and the Poisson 
bracket implies that the two brackets combined define a Lie-Jordan algebra. 

In particular, by considering generic quadratic functions of two complex coordi- 
nates, we find a complex valued quadratic function whose real and imaginary parts 
are quadratic functions of the previous type. All in all, the result is: 

Lemma 5. Complex valued quadratic functions close on a C* -algebra with respect 
to the Hermitian bracket. 

Thus we have found that Hermitian operators are associated with Hamiltonian 
vector fields which are also Killing. As a matter of fact, this property characterizes 
functions on C 2 which are associated with Hermitian operators. 

5.2. The projective space. We can also consider the Hermitian tensor field re- 
lated to the one on the complex projective space. It is important to remark that 
while forms can not be projected, contravariant tensor fields can. This is the reason 
why we introduced the contravariant tensors A and G. Thus by considering 



Because computations on the vector space are easier to carry on with respect 
to computations on the complex projective space, it is convenient to consider ten- 
sors on 7i which reproduce the same results that we would have on the complex 
projective space. Thus we want a tensor Gh defined by 



{/, {<?, h} + } = {{/, g}, h} + + {g, {/, h}}+ 



G = 




we can define a projectable tensor as: 




(38) 



G H (n*df,7r*dg) = G(n*df,ir*dg), 
and identically vanishing on 

(ip\dip) = qjdq j +p dp J + i^dp 3 - p j dq j ). 
The tensor is simple to construct as 



g h = G - (A ® A + r ® r) - i(A <g> r - r <g> A). 



(39) 



Analogously we can introduce 




(40) 
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The next step is to consider the projectable quadratic functions. If we consider 
the basis of the Hermitian operators given by the Pauli matrices, we find: 

_ , fflff2 +P1P2 

^ ql + ql+pl+pl 

QiP2 - Piq2 _ q\+p\-q\- pi 



a2 ql + ql+pl+pl " 3 ql + ql+pl+pl 

We find that only the functions associated with {a\, (72,(73} define non-trivial 

functions on the complex projective space. Of course, their associated vector fields 

generate the algebra of SU(2). 

We can compute now the action of the tensor G on these functions, and we 

obtain: 

G(e CTo ,/)=0 V/ 
G(e CTl ,e CTl ) = e - 4e^ 

G(e CT2 ,e CT2 ) = e - 4e^ 2 

G(e a3 ,e a3 ) = 4(e - e\) 

<5(e CTl ,e CT2 ) = -4(e ai e a2 ) 

In an analogous way, other products can be computed. But we already see that, 
due to the denominator, we do not find "orthogonality" of the expectation values 
associated with different Pauli matrices. Instead, what we obtain is 

Lemma 6. The action of G on the set of projectable functions corresponds to 

G(e A ,e B ) = e [A , B]+ - e A .e B . 
This implies that for A = B we have 

G(e A ,e B ) = e A 2 - e\, 
i.e. we find the variance, the quadratic deviation from the mean value. 

As a conclusion we get 
Corollary 1. G is directly related to the indetermination relations. 

6. Example II: three level quantum systems 

6.1. The projective space. In the case of a three level quantum system, the 
Hilbert space description is done on C 3 , the operators being (up to multiplication 
by the imaginary unit) elements of the Lie algebra u(3). 

The description at the Hilbert space level is fairly similar to the one corre- 
sponding to the two levels system. Now the orthonormal basis we consider is 
{| e i)> I e 2 } , |e 3 )} and the corresponding coordinates, complex and real will be {zi,z 2 , z 3 } 
and {qi,q2,q3,Pi,P2,P3}- The two restrictions: 

= 1 IV) e c 3 

and 

ipi ~ V2 V2 = AV>i A e C, 
define now the corresponding sphere on the projective space VC 
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We can read from these relations the definition of the complex and real coordi- 
nates and the subsequent expressions for the tensors. But instead of repeating this 
line of argument, we shall use a different approach and use the momentum mapping 
H : H — > u*(3) to describe the system. 

As we saw above, the relevant structures now shall be those of u*(3), which are 
defined as: 

• the Poisson structure corresponds to the canonical Lie-Poisson structure of 
the dual of a Lie algebra, defined as: 

{a u , a v } = a [U)t) ] Vu, v e u(3), (41) 

where we denote by a u (resp. a v ) the clement of u*(3) corresponding to 
the element u e u(3) (resp. v <G u(3)). 

• the Riemannian structure is associated to the canonical Jordan structure 
of the Lie algebra u(3). Indeed, we define the tensor R 

TZ(0(da u , da v ) = £([«, v]+) = l -T^{uv + vu) , (42) 

Thus we see that in order to compute the expression of these tensors in some 
basis, we have to compute the structure constants of the u(3) algebra. Thus if we 
take as basis the one defined by the Gell-Mann matrices 



( 


I— 1 






Ai = 1 





i) 


A 2 = | 


V o 









(° 









A 4 = 





:) 


A 5 = 


I— 1 












A 7 = | -i I A 8 = — | 1 | A = \/-| 1 

10 1 

They satisfy the scalar product relations 

TrA^A„ = 28^ v 

Their commutation and anti-commutation relations are written in terms of the 
antisymmetric structure constants and symmetric d-symbols d M „ p . We find 

[Ap, A^] = 2iC^ vp \p [Ap, A p ] + = 2W -Ao^^ + 2d iivp \ p . 

The numerical values turn out to be 

\/3 1 

Cl23 = 1) C458 = C678 = ~^-Cl47 = ~ Cl56 = ^246 = C257 = C345 = — C367 = ^ 

The values of these symbols show the different embeddings of SU{2) into 5/7(3) C 
U(3). For the other coefficients we have 
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djjo = -d jj = -djoj = y g 3 = 1 ' ' ' ' ' 8 

— ds88 = djjs = dj S j = j = 1, 2, 3 
d s w - = djjs = d jSj = - j = 4, 5, 6, 7 
rf 3j 3 = 4«3 = dj3j i = I j = 4, 5 d 3j j = = d j3j = ~2 j = 6, 7 

^146 = ^157 = ^164 = ^175 = — ^247 = ^256 = ^265 = ~<^274 = \ 
^416 = ^^427 = ^461 = ^^472 = ^517 = ^526 = ^562 = ^571 = \ 
^614 = ^625 = ^641 = ^652 = ^715 = ~ ^724 = ^751 = ~ ^742 = \ 

Thus we have completely characterized the Ricmannian and the Poisson tensors 
along with their associated brackets on the density states of the three-levels system. 

7. Example III: a bipartite composite system of two spin \ particles 

Our last example will be devoted to the study of a system composed of two spin- 
1/2 particles. In this case the Hilbert space is the 4 dimensional complex space 
C 4 = C 2 ® C 2 and the set of Hermitian operators corresponds to the elements of 
the Lie algebra u(4) after multiplication by the imaginary unit. 

The importance of the example arises from the fact that this is the simplest 
situation where entanglement of two systems can be studied. 

We will use then the approach used in the previous section, and describe the 
geometrical structure of the composite system by using the momentum map and 
thus the image of the Hilbert space C 4 on u* (4) . 

7.1. The choice of coordinates on u*(4). The most suitable choices corresponds 
to the tensor product of the basis of the two u(2) factors corresponding to each 
subsystem: 

B={o , fc®7j-}fcj=0,l,2,3' ( 43 ) 

where o-q stands for the identity matrix in dimension 2 and {ai}i=i t 2,a and {71)1=1,2,3 
represent two copies of the Pauli matrices. Thus we represent the elements in u(4) 
by a set of 16 real coordinates {A , m k , nj,r k j} for k, j = 1, 2, 3 defined as: 

333 
u(4) 3 A = \o^o®To + ^in k cr k (g)T a + ^n k (7 a (g)T k + ^ ('n/s%+r fej )a-/ s ®rj (44) 

fe=l fe=l k,j = l 

Now all tensor fields shall be written in terms of these coordinates in the next 
sections. This basis has an immediate extension to u*(4) by the duality relation 
defined by the invertibility of the Killing-Cartan form. The main advantage of this 
choice of coordinates is that separable elements in the space of density operators 
correspond to the points defined by setting 

r o =0 Vi,j (45) 

7.2. The Riemann and the Poisson structures. The computation of the Pois- 
son and the Riemann tensor requires again to compute the structure constants of 
the Lie and the Jordan algebras in this basis. Let us present the results by means 
of the corresponding Hamiltonian and Ricmannian vector fields associated to the 
basis elements: 
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7.2.1. Hamiltonian vector fields. Besides the Hamiltonian vector field associated 
to A , which is a Casimir of the structure, these are the Hamiltonian vector fields 
associated to the elements of the basis. We list just a few of them, since the others 
can be obtained from these ones by a simple permutation of indices 



{mi,-} = r 33 - r 23 - 1- ^32 ^ r 22 - h 

dr 23 dr 33 dr 22 dr 32 

9 9 9 9 

r-21 r 3 i ^ m 2 h m 3 - — 

arzi or 2 \ am 3 drri2 

, , d d d d 

{ni,-\ = -f327 l" r 33^ r 22 - h r-23 ^ 

or 33 dr 32 dr 23 ar 22 

d d d d 

ri2^ h ri3- n 2 - h n 3 - — 

ori 3 or\2 dn 3 dn 2 

{rn,-} = ri2Tr + m D n ? + 2miri 2 ) tt 

dn 3 vv ' ' dr 13 

d d 

ri 3 T{ 1- ((-1 + ml) n 3 + 2m 1 r 13 ) h 

dn 2 vv ' ' dr 12 

9 it-, 2\ \ & 

r 2 i^r (m 2 (-1 + nj + 2mr 2 i) 

dm 3 dr 3 \ 

9 ii, 2 \ \ 9 

r 3 \ ~ 1- ( m 3 { - 1 + ri! j + 2nir 3 i j h 

am 2 c*r 2 i 

(m 2 ri 3 + mi (m 2 n 3 + r 23 ) + n 2 (m 3 m + r 3 i) + nir 32 ) 

dr 22 

d 

(m 3 r 12 + n 3 (m 2 ni + r 2 i) + mr 23 + mi (m 3 n 2 + r 32 )) h 

Cr 33 

(-m 2 ri2 - mi (m 2 n 2 + r 22 ) + n 3 (m 3 m + r 3 i) + mr 33 ) 

or 23 
9 

+ (m 3 n 3 - n 2 (m 2 ni + r 2 i) - nir 22 + mi (m 3 n 3 + r 33 )) - — 

9r 32 

9 9 

{ri2,-} = ((-1+m?) ni + 2mini) 

dn 3 dr 13 

9 ii ., 2 \ „ \ 9 

n 3 ^ h ( (-1 + mij n 3 + 2miri 3 j 

oni vv ' ' arn 

-r 22 t^- + (m 2 (-1 + ra|) + 2n 2 r 22 ) + 

, i \ i w 9 9 

(-m 3 ru + n 3 (m 2 n 2 + r 22 ) + n 2 r 23 - mi (m 3 m + r 3 i)J h r 32 - h 

or 33 om 2 

(m 2 r i3 + mi (m 2 n 3 + r 23 ) - n 2 (m 3 m + r 3 i) - mr 32 ) 

(m 3 (-1 + n 2 ) + 2n 2 r 32 ) 

9 

(m 2 rn + mi (m 2 m + r 2 i) + n 3 (m 3 n 2 + r 32 ) + n 2 r 33 ) h 

or 23 

(m 3 ri 3 + n 2 (m 2 ni + r 2 i) + njr 22 + mi (rn 3 n 3 + r 33 )) 

or 3 i 
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d d 
{r2i,-} = r n^> ( m i (-1 + n i) + 2nim) 

7*227^ + ((-1 + ml) 7i2 + 2m 2 r 2 2) 77^ h 

7*23 TT^ ((-1 + Til) 713 + 2m 2 r2 3 ) 77^ 

0712 07*22 

9 / / , 2\ \ 9 
7-3177 h (?Ti3 (-1 + n 1 ) + 2nir 3 ij h 

07711 07'n 

(-m 2 ri3 - mi (m 2 n 3 + 7-23) + 7i 2 (m 3 ni + r 3 i) + nir 32 ) 77 h 

Cri2 

d 

(-n 3 (mini + m) - mri 3 + ra 3 r 22 + m 2 (m 3 n 2 + r 32 )) 77^- + 

(m 2 ri2 + mi (m 2 7i2 + r 2 2) + n 3 (m 3 m + r 3 i) + nir 33 ) 

0ri 3 



(n 2 (mim + m) + mri 2 + m 3 r 23 + m 2 (m 3 n 3 + r 33 )) - — 

07"32 



7.2.2. Riemannian vector fields. Also by using the Riemannian structure we can associate 
vector fields with coordinate functions. Again, we just list a few of them: 



9 9 9 9 9 d d 

J\ a = 711- 1-712 77 1~ 71 3 77 1-771177 h 7712 77 h m 3 77 h 7/0 77 h 

0ni 0n2 0n 3 dmi 9m,2 0013 oyo 

d d d d d d 

7-11 77 h 7-12 77 1" 7-13 77 h 7-21 77 h 7"22 77 h 7"23 77 h 

0rn 0ri 2 0ri 3 0r 2 i 9r 2 2 07-23 

d d d 9 9 9 
rzi 77 h 7-32 77 1- 7-33 77 nmii n 2 mi n 3 mi 7: h 

07-31 Or 3 2 07-33 07-11 OT 12 07*13 

9 9 9 9 9 9 
-nim 2 77 n 2 m2 77 n 3 m 2 77 h -nim 3 - n 2 m 3 - n 3 m 3 77 

07-21 07-22 07-23 Or 31 Or 32 Or 33 

d , . 9 , . 9 , . 9 

4i = 2/0 77 1- (mmi + ni) h (mm 2 + ri 2 ) h mm 3 + m) 77 

ami 9m 9n 2 9n 3 



(mi r-ii + ni (m? + y - l)) 77 (miri 2 + n 2 (m? + y - l)) 77 

v v orn v v " Ori2 

(miri 3 + n 3 (m? + 7/0 - l)) -Jr 

07-13 

1 \ 9 , ,9 

m 2 (mini + 7*11 ) 77 m 2 (mm 2 + n 2 J 77 

9r 2 i OT22 

/• \ , N 9 

m 2 (mm 3 + n 3 ) 77 rn 3 (mmi + mj 77 

or 23 or 31 

, \ 9 , ,9 

m 3 (min 2 + ri 2 ) 77 rn 3 (mm 3 + n 3 J 77 h 

07-32 9r 33 
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/ \ 9 , . 9 

Jni = mmi+m h. h (m 2 m + r 2 i) ~ h 

omi 9m 2 

(m 3 ni + r 3 i) -7, 1- J/o^ (mm + mi (n? + y - l)) 

9m 3 ani v v " orn 

/ n 9 , i 9 

n 2 (mini + r n n 3 (mini + m) 

9n 2 dri3 

(nir 2 i + m 2 (n? + y - l)) "-2 (m 2 m + r 2i ) n 3 (m 2 ni + r 2 i) -= 

v v " or 21 or 22 or 23 

i i i ,w 9 , .9 . . 9 9 

(nir 3 i + m 3 (ni + j/o - 1) ) ^ "-2 (m 3 m + r 3 i) 7; fi3 (m 3 m + r 3 i) -7- h m- — 

9r 3 i 9r 32 Or 33 oy 

Jr-n = (m-irn + m (rn? + y - 1)) + mi (m 2 m + r 2 i) + 

mi (m 3 ni + r 3 i) h (mm + mi (n\ + y Q - l)) h 

9m 3 v v " 9ni 

/ \ 9 , .9 

ni (min 2 + r 12 ) h n x (min 3 + ri 3 ) 

9n 2 9n 3 

((2n? + y - 2) m? + 2mriirai + n? (j/ - 2) + yo) ^ 

(ram 2 r-ii + m (min 2 + n 2 (2m? + y - 2))) - 

9 

(mm 3 rii + rn (min 3 + n 3 (2m? + y - 2))) — 

9 

(ra 2 mrii + mi (rnr 2 i + m 2 (2n x + y - 2))) — h 

(m 3 n 3 — m 2 mn 2 — mm 2 (2m 2 m + r 2 i) + r 33 ) — 

or 22 

9 

(m 3 n 2 + m 2 niri 3 + m x n 3 (2m 2 n x + r 2 i) + r 32 ) 

or 23 

9 

(m 3 nim + mi (mr 3 i + m 3 (2nf + y - 2))) — 

9 

(m 2 n 3 + m 3 mri 2 + r 23 + min 2 (2m 3 ni + r 3 i)) 7; h 

9r 32 

(m 2 n 2 - m 3 mri 3 + r 22 - mm 3 (2m 3 ni + r 3 i)) h (mmi - rn) - — 

9r 33 9y 
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(miri2 + n 2 (m? + y - l)) + mi (m 2 n 2 + r 22 ) + 

mi (m 3 n 2 + r 32 ) h n 2 (mini + rii) h 

am 3 am 

12 ,vv 9 , . d 
n 2 r 12 +mi (n 2 +Vo - 1) ) 1- "2 mi«3 + fi 3 ) 



min 2 r-n + n\ (mma + n 2 (2m? + y — 2))) 



(2n 2 + - 2) m? + 2n 2 n 2 mi + n| (y - 2) + y ) 7^7- 

mimn-i + n 2 (mina + n 3 (2m? + y - 2))) 

d 

m 3 n 3 + m 2 n 2 rn + mini (2m 2 n 2 + r 22 ) + r 33 ) 

or 2i 

m 2 n 2 n 2 + mi (n 2 r 22 + m 2 (2n 2 + 2/0-2))) — h 

d 

m 3 ni — m 2 n 2 r i3 — m x n 3 (2m 2 n 2 + r 22 ) + r 3i ) h 

or 23 



m 2 n 3 - m 3 n 2 rn + r 23 - mini (2m 3 n 2 + r 32 )) 
m 3 n 2 r i2 + mi (n 2 r 32 + m 3 (2n| + j/o - 2))) 



<9r 3 i 
d 
dr 3 2 
d 



d 

m 2 ni + m 3 n 2 r i3 + r 2i + min 3 (2m 3 n 2 + r 32 )) h (mm 2 — r i2 ) - — 

9r 33 ctyo 

By direct computation, but also form general considerations, these vector fields generate 
the Lie algebra of GL(3,C). This algebra contains Hamiltonian vector fields and vector 
fields which describe dissipation or decoherence. These results are related to those of 
Rajeev ([ST]). 

7.3. The analysis of entanglement. Entanglement is a property of composite physical 
systems which plays a very important role in many different phenomena, but in particular, 
it has become a crucial issue of quantum computation and quantum information theory. 
Despite the growing interest in recent years, it was already discussed by Schrodinger and 
the "founding fathers" of quantum theory in the early years of quantum theory f [321 133] ). 
Roughly speaking, entanglement is the situation complementary to separability 

Definition 3. Let be state of a Hilbert space Ti = Ti\ ® TL2 of a bipartite system. 
Then, \ip) is said to be separable if there exists a pair of states |i/>i) £ TL\ and |^/> 2 ) 6 H2 
satisfying \ip) = ® \ip2) ■ 

Thus a state which is not separable is said to be entangled. But entanglement exhibits 
many interesting properties, for instance the fact that there is a gradation in the level 
of entanglement of the different states. Thus we can measure the entanglement of a 
state by using physical quantities. These different observables are called entanglement 
witnesses. 

Usual choices to describe entanglement are the concurrence of the state, the von 
Neumann entropy of one of its partial traces (i.e. the entropy of the density state 
pi = Tr 2 p4, or of p 2 = Trip,/,, where p^ = 

Definition 4. The concurrence of a density matrix p 6 T>(TC) is defined as 

C(p) = max(0, 2\ max (p) - Tr(p)) 

where p corresponds to 

P : 



\J (<T 2 ® (T 2 )p*(02 ® 2 )p 
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and X m ax(p) stands for its largest eigenvalue. 

Definition 5. The von Neumann entropy of a density matrix p G T>(7i) is defined as 

S(p) = Trplog(p). (46) 

When the density matrix corresponds to a pure state, though, the function above vanishes. 
Thus we define the corresponding entropy as the value of the function on the partial trace 
over one of the subsystems: 

S{ Pi> ) = Trpi log(pi) pi = Tn/ty (47) 

On the other hand, if we consider the case of mixed states, the situation is not equally 
simple. In general, it is necessary to consider more than one entanglement witnesses in 
order to completely characterize the state of the system. One interesting question arises 
thus: how can we describe the independence of the different quantities we use? 

In the framework of classical mechanics this question is simple to answer. Given two 
physical quantities, which are represented by two functions /i, f^ on phase-space, they are 
said to be independent at a point p if their exterior differentials satisfy 

(4f 1 A4fa)(p)^0 

The usual approach to Quantum Mechanics, in terms of Hilbert spaces or C*-algebras 
does not allow a similar treatment of the analogous quantum problem. We lack of a 
well developed non-commutative differential calculus which would allow to define a "non- 
commutative" exterior differential and therefore to extend the previous definition to the 
quantum setting [271 1341 [35l |4] 

However the geometrical formalism we introduced in the previous sections allows us 
to look at the problem from a different perspective. Treating the quantum state space as 
a real differential manifold, we do have a differential calculus at our disposal: the usual 
differential calculus of real manifolds. 

Consider the Hilbert space TL and an operator A. We know that we can associate with 
A the quadratic function 

A-*f A (il>) = W\A\i>) i^en. 

In the geometric description of Quantum Mechanics we read from the set of quadratic 
functions the algebraic structures available on the set of operators: 

• the associative product of operators is translated into the non-local product *, 

• the Lie algebra defined by the commutator is translated into the Poisson algebra 
defined by the tensor A 

• the Jordan algebra given by the anticommutator is translated into the Jordan 
algebra defined by the tensor G 

But the geometric description also includes a pointwise algebra (/a-/s)(V') = /a ('!/') /b(V0> 
which is commutative, and whose differential calculus is the standard one. This is the al- 
gebraic structure with respect to which we define the differential algebra we are interested 
in: 

Definition 6. Two observables A and B are said to be functionally independent if 

their associated functions satisfy 

oJa A ojb /0 on a dense submanifold of Ti. 

This is a condition which is simple to test and exhibits a clear advantage offered by the 
geometric approach to the description of quantum systems, for there is no simple analogue 
of this structure in the usual framework. 

Let us see how this works in a particular example where we can test the independence 
of von Neumann entropy and the concurrence of our two qubit system: 
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Example 1. Now we will test this formalism with a particular example. Consider for 
instance the family of density states defined by the matrices: 



Pt = 



as it is considered in |12j . Such a matrix represents a density state provided that 

< a + 6 < 1 < c < 1 Aab>c 

We can take an adapted basis for 
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V o 
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This is clearly a 4- dimensional submanifold S ofD{H). 
it, considering the matrices 
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l/\000 -l/\00 
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We can use the four real numbers {a, b, c, <f>} as adapted coordinates on that submanifold. 

Now we can evaluate the three functions above on these states. As we already know 
the expression of the Jordan and the Poisson bracket we can also obtain the corresponding 
Hamiltonian and gradient vector fields. And besides, we can also study the independence 
of the functions by evaluating the expressions of 

dSAdC 

where d represents the exterior differential associated with the differential manifold con- 
stituted by u* (4) . 

Let us thus proceed: 

• The value of the different functions is easy to obtain. We have that von Neumann 
entropy reads, 

2S{p t ) = -2(-l + a + 6) logfl - a - b] + 



(a + b - y/(a - &) 2 + c 2 ) log i (a + b - ^(a-b) 2 +c^ 
(a + b+ a/ (a - ft) 2 + c 2 ) log i (a + b + ^(a-b) 2 +c 2 



+ 



The value of the concurrence is very simple: 

C(jh) = c (48) 

We can study now the Poisson brackets corresponding to them. It is simple to 
prove that these functions commute, i.e. 

{S,C} = (49) 

This implies that the local transformations generated by them commute. 
Finally, we can study the independence of the different functions. This is an im- 
portant issue, in particular the independence of the von Neumann entropy and the 
concurrence, because it affects the description of entanglement of general quantum 
density states. We can prove the following: 

Lemma 7. The concurrence and the von Neumann entropy of the family of states 
p t are not independent in all the space of density states. 

Proof. We are considering the submanifold of u(Tt) corresponding to the family 
of density states pt- 

Thus the differential of the concurrence is trivial to obtain: 



dC(pt) = dc. 



(50) 
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The computation of the differential of the von Neumann entropy is quite more 
involved. It is evident from the expression above that the functions 5 depends 
on the three variables. But as C depends only on c, we have to consider only the 
a and b dependence in what regards the computation of dS A dC. To this aim we 
compute: 



as as 

da ' db 

Now, the condition for ( |50[ ) to be equal to zero corresponds to 

as _ _ as 

da db 
And these conditions become 

„2 

2Log[l - a - b] + Loi 



ab — 



= Log [4ab - c] - 2Log \a + b + y/(a - b) 2 + c 2 ] = 
These equations have a solution on 

b = a; c = \/-4 + 8a - 4a 2 + 8fe - 4a6 - 46 2 
If we represent the condition for c as a function of a and 6 we verify that it is 



1 1 

3 <a< 2 ; 



well defined for all values of a and b (we represent the function c = ^J— 4 + 8a — 4a 2 + 8b - 
and the function c = y— a + 16a — 12a 2 , which corresponds to the evaluation on 
the submanifold a = 6). 



4a6 - 4o 2 



0J0 




Thus we conclude that there is a nonempty subset of u*(4) where the von 
Neumann entropy functions and the concurrence function introduced above are 
not functionally independent. On any point outside this submanifold the two 
functions are indeed functionally independent. □ 

8. Conclusions and outlook 

To summarize the exposition we would like to mention a few advantages and some 
potential applications of the geometric approach to Quantum Mechanics that we just 
presented to the field of Quantum Information Theory and Quantum Control. 
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The main idea the paper aims to transmit is that Quantum Mechanics can be efficiently 
presented as a geometric theory, considering that a quantum dynamical system can be 
defined as 

• a dynamical system on a Kahler manifold where the dynamical vector field is 
Hamiltonian with respect to the canonical symplectic structure and is a Killing 
vector for the Riemannian metric. 

• a dynamical system defined as a derivation for an associative algebra endowed 
with a pair of operations which define a Lie- Jordan algebra on it. 

In both cases, the message is simple: the geometric framework provides an approach 
to quantum mechanics which is very close to the usual geometric framework of classical 
mechanics and many of the "classical tools" can therefore be used in the quantum domain. 
The most relevant is perhaps the differential structure encoded in the exterior differential 
calculus. Such a structure has not a simple analogue in the usual quantum description, 
but provides a powerful tool to study important notions as the independence of operators, 
as we saw in the last section. Also the passage to the projective space turns out to be 
a simple task for, up to some conformal factors whose meaning is clear in geometrical 
terms, the treatment of the projective theory is analogous to the one defined at the level 
of Hilbert spaces. 

But the analogies of the quantum and the classical frameworks have also some other 
deeper implications, particularly in the realm of Quantum Computation and Quantum 
Information theory. We have seen how the cases of one and two qubits and one qutrit can 
be perfectly described within the geometric approach in fairly simple terms. The geometric 
implications of entanglement is under study now, and we hope to be able to provide a 
detailed description in a short time. But also quantum channels, their capacities, the 
study of decoherence, etc are concepts with a simple translation in the geometric domain. 
We are presently addressing these issues and they will appear elsewhere. 

Finally we would like to mention another advantage of the geometric description of 
quantum mechanics: the study of quantum control problems. Geometric control theory 
is a well stablished theory which has provided very powerful tools to describe classical 
control systems. On the other hand, the usual description of quantum control problems 
makes difficult to translate or adapt those classical tools to the case of quantum systems. 
But once we have proved that quantum dynamical systems can be described using the 
same mathematics which is used in the classical case, the possibilities for application of 
classical geometric control theory seem to be quite rich. And this is even true for systems 
and phenomena purely quantum in origin, as the quantum Zeno effect, whose applications 
to control problems are becoming an interesting topic 38. 1191 fT5] . 
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